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Abstract 

The low temperature physics of disordered systems is governed by the statistics of extremely low 
energy states. It is thus rather important to discuss the possible universality classes for extreme 
value statistics. We compare the usual probabilistic classification to the results of the replica 
approach. We show in detail that one class of independent variables corresponds exactly to the 
so-called one step replica symmetry breaking solution in the replica language. This universality 
class holds if the correlations are sufficiently weak. We discuss the relation between the statistics 
of extremes and the problem of Burgers turbulence in decay. 
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1 Introduction 

The rephca method is one of the very few general analytical methods to investigate disordered systems 
[p. Although the physical meaning of Parisi's 'replica symmetry breaking' (rsb) scheme needed to 
obtain the correct low temperature solution of various random models has already been discussed on 
several occasions [Q] , its precise relation with the so-called extreme value statistics [||, ^ (and therefore 
its scope and limitations) was not previously clearly established. That such a relation should exist is 
however intuitively obvious: at low temperatures, a disordered system will preferentially occupy its 
low energy states, which are random variables because of the disordered nature of the problem. The 
statistics of the free-energy (or of other observables, such as energy harriers Q) will thus reflect the 
statistics of these low energy (extreme) states. It is well known in probability theory that extreme 
value statistics can be classified into different universality classes [Q, ||. Conversely, the RSB scheme 
has shown the existence of at least two broad classes of systems, those with a first order, 'one step' 
RSB and those with continuous RSB. 

It is easy to identify the 'one step RSb' class with one particular universality class of extreme 
value statistics, i.e. the Gumbel class, which concerns the minimum of continuous variables which are 
unbounded but have a distribution decaying faster than any power at —00. The simplest representative 

* Unite propre du CNRS, associee a I'Ecole Normale Superieure et a TUniversite de Paris Sud 



1 



of this class is the Random Energy Model (rem) |^], where the energy states are independent (but 
not necessarily Gaussian distributed). An interesting point is that the REM can be given a spatial 
structure, for which the replica theory still provides the exact solution. This spatial rem is in turn 
connected, in one dimension, to the problem of decaying Burgers' turbulence in which an infinitely 
compressible fluid evolves from random initial conditions. Exact results for the velocity correlations 
at large times in Burgers' turbulence have been obtained long ago by Kida [^. We shall show that 
these results coincide with those of the replica method, the underlying reason being that the late stage 
of turbulence decay is governed by the extreme values for the integral of the initial velocity field. 

It is less easy to identify the other universality classes of extreme value statistics. There should be 
at least two types of generalizations. One type still concerns independent random variables but with 
either power law decay of the distribution (in which case there is a priori no replica formalism), or 
bounded random variables (the Weibull distribution of extremes), which does not seem to correspond 
to any known RSB scheme. The other type concerns correlated variables, for which the only results 
known to us are actually derived in the framework of replicas: those are cases of full RSB, which 
describe random variables with a certain (hierarchical) type of correlations. 

These universality classes are the counterpart for extremes of random variables to the usual univer- 
sality classes studied in the framework of sums of random variables. Taking the well known example 
of random walks or polymers, the usual random walk, or ideal polymer, is described asymptotically by 
the Gaussian central limit theorem, while the addition of independent variables with infinite variance 
leads to new universality classes (Levy sums) [|^]. The introduction of long-range correlations like 
self avoidance also leads to totally new universality classes fl^. We wish here to take a first step in 
an analogous categorization for extreme values, which appear naturally in disordered systems at low 
temperatures. 



2 Extreme value statistics 



2.1 Scaling regime 

We start by recalling standard results of extreme values statistics, in order to set the stage for the 
following discussions. Consider M independent, identically distributed random variables Ei, i = 
1,...,M ('energies'), such that the probability distribution decays for Ei — > — oo faster than any 
power-law: 

P{E) -^exp[-B\E\^] B,5>0; E^-oo (1) 

We are interested in the statistics of the lowest energy state E* = minjii^i, ...Em} for large M. Defining 
'P< (E) as the repartition function of E: 

rE 

V<{E) = / dE'P{E') (2) 
one can express the distribution Pm of E* as: 

Pm{E*) = MP{E*)[1-V^{E*)]^'-' = -^[Vy{E*)]^ (3) 

For large M, the minimum E* will be negative and large, so that: 

[1 - V<{E*)]^^ ~ exp[-MP<(^*)] (4) 
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The repartition function of E* thus becomes very small when E* is smaller than the characteristic 
value of the energy Ec defined by MV<{Ec) = 1. To logarithmic accuracy, this gives in the case of 
the distribution (|l]): 



B 

Defining now E* = Ec + e, with e <C \Ec\, one has, to first order: 

[1 - V^{E*)]'' exp[-exp(S5|S,|^-^e)] (6) 

Finally, introducing the rescaled variable u = B6\Ec\^~^ e, one finds that the rescaled variable obeys, 
for large M, a universal 'Gumbel' distribution ^ P*, independent of the coefficients A,B and 
exponents a, 6: 

P*{u) = exp(u — exp-u) (7) 

A very important property, which we shall emphasize later on, is that P*{u) vanishes exponentially 
for u — > — oo (and much faster still for u — > +oo). The maximum of P*{u) occurs at u = 0, meaning 
that Ec is actually the most probable value for the extreme energy. Finally, as in any 'central' limit 
theorem, this behaviour is only valid in the region where the deviation e from Ec is of the order of 
E^-^/B, which goes to zero with M if 5 > 1 and diverges otherwise. The relative fluctuations e/Ec, 
however, are always of order 1/logM. 

2.2 The large M limit and the Random Energy Model 

Let us now consider the following partition function: 

^ E- 
Z = ^Zi Zi = exp[--^] (8) 

1=1 

where the Ei are distributed as in (|^). This is a slight generalization of Derrida's original rem, initially 
introduced with a purely Gaussian distribution (6 = 2). Obviously, the independent variables Zj, are 
large when Ei is large and negative. In the scaling region defined above, due to the exponential tail 
of d^, the distribution of z decays for large z as a power-law: 

P{z) oc z-^-f" ■ z^oo (9) 

where Ec is the most probable ground state energy of the system, given by Eq. (|5|), and 

/i = TB6\Ec\^-^ . (10) 

The partition sum Z behaves very differently in the region fj, < 1, where the average value of z diverges 
and thus only a small number of terms (those of order M^^^) contribute to Z, and in the region /U > 1, 
where all the M terms give a (small) contribution to Z. This means that for 

Tc = (11) 

for which fi = 1, the probability measure concentrates onto a finite number of states, corresponding to 
the glass transition in these models. In the random energy model, M is the number of states M = 2^ . 
In order to have an extensive ground state energy (Ec oc N) and Tc finite in the large N limit, one 
should choose (see Eqs. (|^JTl|)) B = N^~^. For S = 2, this indeed coincides with the usual scaling of 
the energies in the rem. 
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Let us now study the statistics of the weights pi = Zi/Z in the glassy region T < Tc- Since: 

where Z' = J2k(yti) is independent of Zj (and of order M^^^), one readily finds that 

^(-) = TT^P ( - = ^) (13) 



(1 — w)^ \ 1 — w 

For Wi to be non zero in the large M limit, Zi has to be large. In that region one can use the asymptotic 
form (^) for P{z), giving: 

P{w) = -^(1 - w)f''^w-^-^' w > M-V/^ (14) 

where C is a constant fixed by the condition M dwwP{w) = 1. From this probability distribution of 
each weight, one can deduce the moments 1^ = J2i ) which characterize to what extent the measure 
concentrates onto a few states: if all weights are of the same order of magnitude, then ~ M^~'' — > 
for A: > 1; while if only a finite number of weights contribute, the moments Yfc remain finite when 
M — > oo. In the present case, one finds, for /x < 1, 

(see also [jll|). Since fi = T/Tc, one finds that Y2 goes linearly to zero for T ^ Tc, and that 

Yk = i- {r'[k] - r[i])/r[A:])r/r, for r ^ o. 

Finally, the average energy per degree of freedom of the system is constant throughout the low 
temperature phase (T < Tc) and given by 

E/N = Ec/N+ <u> B5Ei-^/N ~ -(log2)i/^ + 0{l/N) , (16) 

where the average < li > is taken over the Gumbel distribution, giving: < u >= r'[l]. 

3 The replica approach 
3.1 The REM 

We shall now show how all these results can be recovered using the replica method. We suppose that 
5 > \ (the case 5 < \ will be discussed below) and introduce the characteristic function g{\) through: 

/oo 
dEP{E) exp[-AS] = exp[5(A)] (17) 
-00 

Since B = N^^^, this integral can be computed at large N with a saddle-point method, which gives: 



g{X) = {5-l)N(-] (18) 



^ We denote as P(.) the probability density of the variable appearing in the parenthesis; hopefully there is not ambiguity 
in the following. 
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In the replica method we need to compute the moments of the Z distribution: 



— Zi^Zi^...Zi^ — exp 

11,12, ...in i\,i2,...in 



1 " 

i a=l 



(19) 



The averaging over the Ei gives: 

Zr^= ^MT.9{7fT.k^a)] (20) 

11,12,. ..in i a=l 

The point now is to understand which configurations of \i\^i2, ■•■in\ will dominate the above sum 
when N ^ oo (and n — > 0). The simplest Ansatz, corresponding to the largest phase space volume, 
assumes that all ia are different, leading to: 

= M{M - 1)...(M -n + 1) exp[n5(^)] ~ exp[n(log M + g{^))] (21) 

Taking n — > 0, one thus finds that the free energy per degree of freedom / = — -^log Z takes the value: 

/ = /o = -riog2-(5-i)5-^r-^ (22) 

The entropy so = —dfo/dT is therefore equal to: 

so = log(2) - (5r)-^ (23) 
and becomes negative below a critical temperature 

T, = ^log(2)^ (24) 

So this solution, called 'replica symmetric' (since all replicas ia play a symmetric role), has to be 
modified in the low temperature phase. The correct configurations which dominate the sum ( [20| ) 
at T < Tc are called 'one step replica symmetry breaking' and are such that the n replica indices 
{ii,i2, ...in} are grouped into n/m groups of m equal indices, which can be written after a proper 
relabelling: 

11 = 12 = ... = im = ki (25) 

im+l = ini+2 = ... = i2m = ^2 (26) 

(27) 

in—m+l — ••• — in — ^n/m (2^) 

and now the indices fci, £^re all different one from the other. These configurations contribute 
to as: 

Z- = M{M - 1)...(M - n/m + 1) exp[-5(-)]— ^ (29) 
from which one immediately deduces: 

/(T) = /o(r/m) (30) 
where /o is defined in Eq. (p2|). The extremum of this free energy with respect to m is obtained when 

df 

^ = = so{T/m) (31) 
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which gives 

T 

m = — = fi (32) 

Note that this relation is independent of 5. Therefore this one step RSB solution predicts that the 
system freezes at the critical temperature Tc which is the temperature where the entropy sq vanishes. 
The energy density is constant throughout the low temperature phase, and equals: 

e = /o(r,) = -(log 2)'/' (33) 

in agreement with the direct computation (|l^). Since the free-energy is constant, the entropy of the 
whole low temperature phase is zero [^]. 

It turns out that also the finer details, like the distribution of the weights of the configurations 
which dominate the low temperature measure, can be computed by this replica approach [p!^ ]. By 
definition, the moments are equal to: 



k 



= ^™n~^ Tw rVT^ J2 Zi^-ZinYlK.Aa, (35) 

\ ' \ I ai,..,ak ii,...,tn j=l 

where the sum primed over the a's runs from 1 to n, with all a's different. Owing to the structure 
of the RSB, this means that one simply has to pick the k < m replica indices ai, ...,ak in the same 
'group', for which there are (m — l)...(m — k + 1) possibilities once oi has been chosen. Hence: 

n{m — l)...{m — k + 1)^—^ T[k — n] 
' = n{n-l)...{n-k + l) ^ r[A;]r[l-M] ^^^^ 



in agreement with the direct computation (|1 
3.2 The REM with S < 1: a first order transition 

The above method fails when 5 < 1, which actually corresponds to a different universality class from 
the point of view of critical phenomena, while the nature of the low temperature phase leaves it in the 
same class as the systems with 5 > 1, in agreement with the extreme value classification which does 
not distinguish between 5 > 1 or 6 < 1. In order to study the transition, we use Derrida's original 
'microcanonical' method. Using the normalisation B = N'^~^ , the partition function is equal to: 

, lei 

Z= deexpNifie) 99(e) = log 2 - lel"^ + ^ (37) 

where e = E/N, and Cc is the energy density beyond which there are no states (for N — > 00), i.e.: 
2^exp(— A^e^) = 1. As shown in Fig. (||), the integral is dominated either by e = or by e = — Cc, 
depending on the temperature. When T > Tc = (log 2)^"'^/'', the free energy is equal to — A^nog2, 
while for T < Tc, the free energy is equal to a constant —Ncc = —N{log2)^^^. The transition at Tc is 
now a first order transition from the thermodynamic point of view, with a jump in the entropy. This 
is in contrast with the usual case 5 > 1 where the transition is thermodynamically of second order 



Although there is a jump in the Edwards- Anderson order parameter 
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Figure 1: </?(e), as defined in Eq. (|37|), as a function of e for different temperatures. The saddle point 
is tlius at e = for T > Tc and at e = — jed for T < T^. Tliere is no states (in the hmit ^ oo 
beyond — |ec|. 
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In the low temperature phase, only the neighbourhood of Ec = —N\ec\ is of importance, and we 
get back to the universal Gumbel distribution since the density of states is still locally exponential: 

P{E = i?, + e) oc exp ^ ^i^-5 (38) 

The value of /i again determines the statistics of the weights, as above. Note that however that for 
(5 < 1, the value of the parameter /i (corresponding to the RSB parameter) is smaller than 1 at the 
transition T = . Hence I2 is discontinuous at T = Tc, in contrast to the case 5 > 1. 



3.3 Physical interpretation of the rephca solution 

The reason why replica symmetry must be broken in order to get sensible results in this problem is 
rather clear. Since the distribution of the Boltzmann weights Zi is a power law with an exponent ^ < 1 
in the low temperature phase, all integer moments of Z (and thus of z) are formally divergent, and are 
thus dominated by a cut-off for large z which has nothing to do with the value of the 'typical' z's, and 
hence of the free-energy. Calculation based on a simple analytic continuation of the results obtained 
for n > 1 are thus bound to fail. The replica method with one step RSB manages to compute {z"^) with 
m = /i, which precisely picks up the contribution of the typical region of z. (Smaller values of m would 
be mostly sensitive to very small z, while larger m's probe atypically large values of z.) The algebra 
corresponding to one step RSB exactly reproduces the extreme value statistics in the case of fastly 
decaying distributions. In this respect, RSB does not mean more than a 'localisation' of weights onto a 
small subset of all configurations, in the sense that major contribution to the partition function comes 
from a finite number of configurations (i.e. all > 0) Actually, the quantities 1^ were also 

introduced in the context of electron localisation in disordered potentials, and called 'participation 
ratios' |13]. 



4 A d-dimensional Random Energy Model 

In this section we want to study a generalized version of the rem, where the energy levels are embedded 
in a euclidean space. Besides its intrinsic interest as a model for a particle in a disordered environment, 
this problem turns out to be also directly relevant to the study of declining Burgers turbulence, as we 
shall discuss in detail in the next section. 

The model is defined as follows. To each point x of a (discretized) d-dimensional space, one assigns 
a potential energy E{x) which is a random number picked up independently on each point, from a 
distribution P{E) the tail of which is given by (1). The total energy on this point is the sum of a 
deterministic part, which we take for instance equal to nx^ /2, and this random contribution E{x). 
This defines a certain energy landscape, to which we associate a partition function Z as: 

Z = j d'^xexp , V{x) ^ ^ + E{x) (39) 

Here we adopt a continuum notation but an ultraviolet cutoff (lattice spacing) is implicitely assumed 
when necessary. The role of the deterministic part proportional to k is twofold. First of all it allows 
one to define a topology in the space of the points x (The limit k = Ocoincides with the rem, the fact 
that the points sit in a d-dimensional space being irrelevant). For this purpose the deterministic part 
could be rather arbitrary, and indeed one can solve the problem with a more general deterministic 
energy. As we shall need a quadratic term later on, and in order to keep the presentation simple, we 
restrict to this particular case. Second, the presence of this confining term allows one to deal with 
this model without the need of introducing a finite box. 
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This model with d = 1 was in fact introduced and studied long ago as a toy model of an interface 



in a random medium |14|: one possible interpretation is that x is the coordinate of a particular point 
on the interface, which feels a random pinning potential E(x), while the quadratic potential is a 
mean-field description of the elasticity due to the rest of the interface. Another interpretation (in the 
context of Bloch walls) is to neglect the deformation of the interface, which is only described by its 
center of mass coordinate x. The quadratic potential is then induced by the demagnetizing fields due 
to the surrounding Bloch walls. 

We want to compute the low temperature properties of this system in the limit when k — > 0. 
For instance one would like to know the typical displacement of the ground state, measured through 
< x'^ >, or the average ground state energy, etc. In the special case where the energy is Gaussian 



distributed, this problem has already been studied by scaling arguments [14|, or with a Gaussian 
replica variational method We shall provide hereafter the exact solution, first using a direct 

extreme value statistics approach and then with the replica method. 

4.1 Extreme value approach 

For simplicity, we restrict to the case d = 1; the extension to higher dimensions is however immediate. 
For temperatures going to zero, we want to find the minimum of all the energies ^x'^ + E[x) when x 
scans a one dimensional lattice. The joint probability that this minimum is achieved on a point x* 
and takes a value V{x) = ^x*'^ + E is given by: 

P{x\E)=P{E) n (l-^<(^ + r*'-r")) (40) 

For K ^ we can safely take a continuum limit and we get: 

P{x\E) = Y^^exp (I dx'log [l - P< [e+^x*' - ^x'^])) (41) 

Integrating over E we get the probability that the minimum is achieved on point x*. For small k, the 
minimum E is expected to be negative and large, and hence only the region where P< is small will be 
of importance. Rescaling x* as x* = x*/y^, we obtain 



P{x*) ~ / dEP{E) exp - / E + (42) 




For small k it is thus clear that the relevant energy region is the one around the value Ec such that 
P<{Ec) = ^, or: 

(Notice that the role of the number M of energy levels in the first section is played here by the length 
scale I/^/kj which is natural.) Expanding the energy around Ec as E = Ec — x*'^/2 + e, we get: 



< 



^^x*^_^\ ^ ^^^p UB\Ec\'-He-^)\ (44) 



The integral over z in (^3|) is a thus a Gaussian integral. We finally get, after a simple integration 
over e: 

P{x*) oc exp l-6B\Ec\^-^^] (45) 
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Therefore we have shown that the typical distance to the origin of the point x* corresponding to a 
minimum energy is 

1 -X 



and more precisely the distribution of x*/^ is a Gaussian of unit variance 0. 

We can also compute the probability distribution of the ground state energy V* as 



(46) 



P{V*) = I dx I dE 6 (V* - -x' - E] P{x, E) 



(47) 



where P{x,E) is given in (41). The result is the following: introducing the rescaled energy u as: 



V* = Ec + 



1 



2Bd\Ec\^ 



B5\E, 



\5-l 



27r 



+ 



B6\E,\^-^ 



(48) 



one finds that u is distributed according to the universal Gumbel distribution, Eq. (0). In particular, 
the extremely deep states are exponentially distributed, as exp[//V*/T], with n = TB6\E^ 



\5-l 



4.2 Replica approach 

Interestingly, the replica approach with a one step RSB also leads to the exact result. Introducing 
again the generating function g{X) of P{E), we have, for large A: 



5(A) 



5 _ 1-1/5-1 /;^x 5/5-1 



B 



(49) 



The average of can thus be expressed as: 

-Z" = ^ exp 



K ^ 1 " 



(50) 



a=l X a=l 

Let US make the ansatz that at low temperature, the saddle point of this expression is such that: 



Xi = X2 
Xm+l = Xm+2 = 



X2m = 2/2 



Xn—m+l — ••• — Xn — Un/m 

and perform the Gaussian integration over the m. We finally obtain: 



Z"' = exp 



n 



m 



1 , 2ttT .m. 
-log +£/( — ) 

2 Km I 



n 



exp--/(p) 



(51) 
(52) 
(53) 
(54) 



(55) 



with p = m/T. Looking for the extremum of / as a function of p we find, in the limit k ^ (and 
with 5 > 1), 

p* = 5^1/5 (log = B5\Ec\^-^ (56) 

where Ec is given by Eq. (I 



^Notice that for small k we have << \Ec\ which justifies our expansion around in the derivation of P{i 
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It is easy to check that the free energy f{p*) precisely reproduces the above result for the ground 
state energy obtained directly, Eq. (|8|). As explained above, the calculation of the quantities Yfc 
within the replica method indicates that the low energy states are exponentially distributed with a 
parameter given by p* = m/T. Hence, comparing (|4^ ) and (|56|), we see that the replica method indeed 
predicts the correct statistics of deep states. The replica method also allows one to calculate 



P{x) = ^ exp 



X2,..X„ 



a=l X a=l 



(57) 



X\=X 



Within the above one step solution, this immediately leads to the following Gaussian result: 



PM=/^exp-^ (58) 

which is identical to Eq. (|45|). 

The replica method also allows one to discuss the non zero temperature regime, which is much 
harder to study directly. As shown above, there is a phase transition towards a 'delocalised' phase 
where = when p = p*T = 1 However, for any non-zero temperature and for 5 > 1, the system 
eventually reaches p = 1 for small enough k. This can be interpreted as follows: as k ^ 0, the number 
of accessible states diverges. But since the difference between the ground state and the first excited 
state decreases as |£'c|^~'' (when 6 > 1), it does become smaller than T for a sufficiently small k, 
beyond which a large number of quasi-degenerate states contribute to the partition function, as in the 
high temperature phase. Only for 5 = 1 is there a true transition temperature, independent of k (see 
[^] for a discussion of this point in a different context). For 5 < 1, one expects a first order phase 
transition (see above). 

Finally, let us note that in the case where the confining potential is harmonic (i.e. equal to kx'^ /2), 
the Gaussian variational replica method developed in ||^, 17] also gives the exact result for p* . 



4.3 Physical interpretation of the rephca solution 



Using the replica method, one can also compute higher moments of P{x), such as P{x)P{y), etc. One 
can then show that the replica solution is identical to the following probabilistic construction for P{x) 
for a given sample: 

P(x) = ^J2^-Sx,x^ (59) 

where the Wa are random weights, chosen with a probability distribution given by Eq. (|l^, and the 
Xa are random variables, independent from the w^s, and chosen according to a Gaussian of width 



5 The random energy model and Burgers' turbulence 
5.1 The Cole-Hopf transformation 

It is well known that the solution of Burgers equation with a random initial velocity field can be 
expressed as a partition sum of the form Eq. (^). Let us restrict for simplicity to one dimension, 
although, again, generalisation to higher dimensions is possible. The Burgers equation in the absence 
of forcing reads: 

dv ^ dv _ d'^v 

dt dx dx"^ 

*Note however that there is a true phase transition only in the hmit «: — > or d ^ oo, i.e. when the number of 
degrees of freedom is infinite. Otherwise, the transition for fj, = 1 is really a crossover. 
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where v is the viscosity. The initial velocity field v{x,t = 0) will be chosen as v{x,t = 0) = ^^j-^- 
Writing v = — 2i/ '^ g^^ , allows to transform the Burgers equation into the following linear diffusion 
equation (Cole-Hopf transform): 

with initial condition Z{x,t = 0) = exp[—E{x)/2i']. The solution thus reads: 

f^°° dx ( I Ax-xof ^, ,\ 



which is, up to a multiplicative factor, identical to the 'spatial' rem defined by Eq.(p9|) with the 
following identification: 

Physically, the disordered problem associated to this rem is that of a point particle interacting with a 
(random) pinning potential E{x), attached by a spring to point xq, which is a simplified model for an 
extended elastic object in a random potential. This model was also recently considered in the context 
of solid friction [p^ j. 

Although the two problems, spatial rem on one hand and decaying Burgers turbulence on the 
other hand, are formally identical, they may differ by the type of questions one wants to address. 
For instance in turbulence one is interested in the correlations of velocities, which involves knowing 
the variations of the free energy of the rem (|6^ ) when xq varies. The case where E{x) is random 
with short range correlations correspond to a short range correlated velocity field v{x,t = 0) with a 
'blue' spectrum (i.e. \v{k,t = 0)p oc k'^, where k is the Fourier variable) and the small viscosity (large 
Reynolds) limit corresponds to small temperature in the associated disordered problem. 



5.2 Cusps and shocks 

In the zero viscosity (or zero temperature) limit, the partition function ( |62[ ) can be evaluated by a 
saddle point method. For a fixed xq, one looks for the value of x* such that k{xo — x*)'^/2 + E{x*) is 
minimum. The saddle point construction Q is graphically explained in Fig. (^), for a simple profile 
E{x). For a given xq, one draws as a function of x the parabola V — k{xo — x)^/2 and looks for the 
minimum value of V, called V*{xq), such that this parabola intersects the curve E{x); calling x* the 
intersection point, the saddle point approximation gives: 

Z(xo, t) exp[-^^] ^ exp ("^[^(^o - x*f + E{x*)]) (64) 



For large values of k, the parabola is very sharp, and there is only one 'optimal' intersection point 
X* for each value of xq; to a first approximation, one thus has 2{xq, t) ~ exp[— i?(2;o) /T]. On the other 
hand for very small k, which corresponds to the large time limit of the decaying Burgers turbulence, 
the parabola V — k{xq — x)"^ /2 is extremely flat and the intersection points will be determined by the 
extreme (negative) values of the potential E[x). In this limit, the statistics of the effective potential 
V* - and thus of the velocity field v{x,t) ~ reflects the statistics of the extreme values of E{x), and 
is thus, to a large degree, universal. Generically, the solution x* depends very weakly on xq and the 
effective potential V*(xo) can thus approximatively be written as: 

V*{xo)c^^{xo-x*f + E{x*) (65) 



12 



Figure 2: Graphical solution of the Burgers equation in the limit of small viscosity, in the neigh- 
bourhood of a cusp. The dashed line is the original potential E{x), while the full line corresponds 
to the effective potential V*{x). The curves actually continue beyond the cusp of V*{x), where one 
metastable and one stable saddle point coexist. 
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with a fixed x*, where E{x*) corresponds to a particularly 'deep' minimum x* of the potential E{x). 
This is the generic situation when one varies xq locally; it corresponds to a velocity field which is 
locally linear: 

v{xo) = '^'^J^^^^ = i^ixo - X*) (66) 

(Remember that by definition dV* /dx* = 0). There exist however exceptional values Xg of xq such 
that the first intersection of the parabola and the curve E{x) appears simultaneously at two points 
x\ < X2'- when xq varies from Xg — eioXs + e, the solution x* jumps from x\ to X2- This corresponds 
to a cusp in the minimum value V* as a function of xq (see Fig. (^)). In the language of Burgers' 
turbulence, this is a shock since the velocity v (which is the derivative of V*) is discontinuous at 

Xq = Xs- 



5.3 Decay from an uncorrelated E{x) configuration: Kida's analysis 

Let us now focus onto the case where E{x) is randomly distributed with a short range correlation, 
and the time t is large, corresponding to a very small k. This limit was studied in detail by Kida 
in the context of Burgers' equation (see also [|^). Let us denote by Xa the various values of the 
intersection points between the parabola and the curve E{x) when one varies xq. After a proper 
coarse graining one can totally forget about the correlations of E{x), and thus the {xa} are randomly 
(Poisson) distributed along the j;-axis. If the distribution of E decays as exp — Ul^^l*^, the extreme 
value statistics tells us that the distribution of Ea = E{xa) is of the Gumbel type. The only delicate 
point is to understand what is the effective number of independent variables, M, appearing in this 
distribution. This number depends on k and is determined self-consistently as follows: as xq departs 
from Xa, at some point (because of the quadratic growing term k{xq — Xa)^ /2) will a better saddle 
point X0 be prefered. Since the width of the Gumbel distribution is given by: 

^ (logM)^ (67) 



(see Eq. (^) above), this sets the order of magnitude of the difference between E^ and Ep, which must 
also be, by definition, of the order of ) . Furthermore, taking the correlation length of the 

potential E{x) to be 1, the effective number of independent variables is given by: 

M = \xp- Xa\ ^ ^=^==(logM)^ (68) 

or, to logarithmic accuracy, and using the correspondance k ^ 1/t, 

M oc \/t(logt)^ (69) 

Note that by definition, M is also the typical distance beween two shocks i{t), which is thus seen 
to grow as t^^^ with logarithmic corrections (these corrections disappear for 6 = 1, where the initial 
potential already possesses the universal exponential tail). This is one of the important results of the 
original analysis of Kida. Furthermore, since the local slope of the velocity is k = 1/t, the maximum 
velocity is of order: 

_ m (logt)^ 

t Vi 

which corresponds to a time dependent Reynolds number : 



(70) 



l-(5 



W^ (logt) s ^^^^ 
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which goes to zero (albeit very slowly) when t — > oo for 6 > 1. This is similar to the above remark 
that for any small temperature, the system goes back into its high temperature phase when k — > 0. 

Using this construction, and the full distribution of the Ea, Kida was able to obtain directly the 
large time behaviour of the two point velocity correlation, v{x)v{x + r), which is a universal function 
once the lengths are expressed in terms of the mean distance between two shocks i{t), and the velocities 
in units of Vm^x 1^]- His result is recalled in the appendix. Let us show how one can obtain precisely the 
same results using the replica method, which in fact provides the ful probability distribution function 
of v{x + r) — v{x). 

5.4 The replica analysis 

Let us first note that Eq. (|65| ) can alternatively (in the limit k, T ^ 0) be written as an infinite sum: 

Zixo) = exp[-V*(xo)/r] =^t.„exp[- ^^"^°^~^"^' ], (72) 

a 

where Xa are Poisson distributed with an arbitrary (see below) linear density a. The Wa are inde- 



pendent random variables again chosen according to the distribution (14), with /i given by (see Eq. 

my- 

/i = T6B^^\log M)" V (73) 

That Eq. (|72|) precisely reproduces Kida's construction comes from the fact that, as T ^ 0, the 
distribution of weights becomes so broad that the sum determining Z{xq) becomes entirely dominated 
by a single term, which is the one which maximizes Wa exp[— k(xo — Xa)^/2]. Again, the corresponding 
Xa switches discontinuously as a function of xq, when another value x^ suddenly takes over. This 
construction is independent of the density a, provided that aM ^ 1 (i.e., in the long time limit). 

The crucial point now is that the explicit construction ([72|) actually gives results which are identical 
to those obtained using a replica representation: 



1 " 

Z{xi)Z{x2)...Z{Xn) ='Yexp[- Rn(a),n{b)XaXb] (74) 

71" a, 6=1 



in the limit n — > 0. In the above expression, ir denotes a permutation of the n replica indices, and the 
Rab matrix is a one step RSB matrix Q with elements Rab = r when a and b are in the same diagonal 
block of size m, Rab = when a and b are in different blocks, and Raa = (1 ~ m)r, enforcing the sum 
rule J2b=i Rab = 0. By equivalent results we mean that one can compute quantities like the average 
probability of being in x, 



P{x) = J = j dx2...dXnZ{x)Z{x2)...Z{Xr. 

or the average of the product of the two probabilities to be in x and y, 



(75) 

n=0 



P{x)P{y) = ^j^j^^^ix'ly ^ I dx3...dxnZ{x)Z{y)Z{xs)...Z{xn: 



(76) 

n=0 



or generalizations thereof, by both methods (here the average is taken over the realization of the 
initial velocity profile). It has been shown in [pO| (and we recall the main steps of the derivation in 



the appendix) that the velocity correlation function v{x, t)v{y, t) can be computed either directly from 
Eq. ([72|), or using the representation (|7^. The important points are the following: 
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- After a proper choice of length and velocity scales, the v{x, t)v{y, t) correlation function is indeed 
identical to the result obtained by Kida (see appendix) , which is the consequence of the fact that both 
approaches actually rely on the universal structure of the extreme events which control the velocity 
field for large times. 

- The present formalism allows us to extend Kida's results in several directions. For example, 
the full probability distribution function of v{x) — v{y) has been computed in |2C]. The problem of 
decaying Burgers turbulence in higher dimensions can also be addressed. 

- The presence of shocks, which manifests itself as a |x — y\ singularity in v{x,t)v{y,t) at short 
distances, is intimately connected with the breaking of replica symmetry |2^, ^ : for a replica symmetric 
matrix Raa = R, Ra^b = Ri, v{x,t)v{y,t) is regular for x — > y. As discussed above, these shocks 
reflect, in the associated disordered problem, the existence of some metastability (see Fig. (|2|). From a 
technical point of view, it is interesting to see on this example how metastability is associated to RSB 
and, as emphasized in to the existence of a short-distance singularity in the effective free-energy 
V*. Precisely the same behaviour is obtained via the Functional Renormalisation Group (frg) |2^: a 
singularity appears in the renormalized correlation function of the effective free energy at scales larger 
than the 'Larkin length', which is the scale beyond which metastability effects become important. 
(However, the way to handle the shocks correctly within the frg is still an open problem ||5|). 



6 Perspectives and other universality classes 

As is the case for the central limit theorem, there are other universality classes, distinct from the 
Gaussian, when one relaxes the hypothesis of a finite variance or of independent variables (or both) 
[p!o|]. This is also true for the statistics of extremes, and it is interesting to discuss how this might 
translate into a replica language. Two main directions can be thought of: independent variables with 
other types of distributions, or correlated variables. 

Let us first consider the case where the energy levels Ei are still independent, but with a tail for 
large negative E decaying as a power-law, j^^l"^"*^. In this case, the extreme values are distributed 
according to the so-called Frechet distribution, which is different from the Gumbel distribution (for 
example, it decays asymptotically as a power-law with the same exponent 6). Rescaling E by M^^^ 
to keep the gap between the ground state and first excited state finite as M — > oo, one can calculate 



the quantities defined in (15). One finds, for M large but finite: 



n = 1 - exp - ( T^j^ ) (<5 < 1) (77) 



,TlogM, 

independently of k. (Similar results are obtained for 6 > 1). This is clearly different form Eq. (|l5|). 
Notice that this case cannot be addressed within the replica method without some modifications since 
all the positive moments Z"-, n > diverge ! 

Another universality class corresponds to Ei which are strictly bounded, i.e. Ei = Eq + e, with 
e > 0. More precisely, the distribution of e for e — > is of the form P{e) = for e small. The resulting 
distribution of extremes is then called the Weibull distribution. Rescaling the energies by a factor 
, we find through a direct computation that is non trivial for all temperatures, i.e., the model 
is always in a low temperature phase. For instance in the case 6 = one gets Yk{T — > oo) ~ kT^~^ , 
and ~ 1 — CT for T — > (C is Euler's constant). This is again clearly different from the 1 step RSB 



result (15). One might hope that such a situation will lead to a new type of RSB, but the situation 



seems more complicated. In the particular case (5 = 0, one finds: 



^"=Eny.^^ (78) 
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where the product is only over the sites i such that X]a=i ^i,ia > 0- The entropy of the rephca 
symmetric solution becomes negative below the temperature Tc = 1/e. Assuming a one step RSB 
saddle point for T < Tc leads to a constant average energy, equal to 1/e; however, the true ground 
state energy can be calculated directly and is equal to 1. Again, this result is not of the replica type. 
In this case however, a sensible replica calculation can be undertaken since all the Z"-, n > are 
convergent. The problem now is that the free-energy cannot be calculated by a saddle point method: 
replica fluctuations are always important. 

In any case, from the point of view of Burgers' turbulence, it is interesting to notice that initial 
conditions for the velocity field which do not belong to the exponential universality class considered by 
Kida will lead to rather different flow structures at long times, even within the class of E(x) functions 
with local correlations. 

Let us now turn to the case where the Ei are Gaussian but long-range correlated; for example the 
case where E depends on a d-dimensional space variable x, and such that: 

mm='-^ (79) 

leading to 

{E{x) - E{y)y (X\X- y|>^ax{0,2-d-,,) ^gQ^j 

The case ij = 0, d = 1 corresponds to a random walk for E{x), which has been studied in detail both 



in the context of Burgers' turbulence |p, and also as a partly solvable spatial rem |1^, 22, The 
general rj, d case has not been solved yet. It has been studied by the Gaussian variational replica 
formalism of Q, which shows [^] that the case rj < 2 — d (corresponding to a growing correlation 
function (^)) requires 'continuous' RSB, while the case rj > 2 — d only requires a 'one step' breaking. 
Independent variables correspond to ry = 2, i.e. a white spectrum for E{q). We conjecture here that 
the case r] > 2 — d belongs to the same (one step RSb) universality class as the rem (rj = 2). It is 
actually not difficult to show directly that the quantities: 

Cn = ^^ n = 2,3,.. (81) 

diverge with the system size below a certain n dependent critical temperature which is independent 
of T] for T] > 2 — d, and identical to those found in the REM. This suggests that the one step solution 
indeed remains exact for all rj > 2 — D. Preliminary numerical simulations ||2^ seem to confirm this. 
This points towards a rather natural result, namely the fact that weak enough correlations (measured 
here by 2 — r/) between the random variables does not change the universality class for the extreme 
value statistics. This is actually a Theorem for the case d = 1: for all r/ > 1, the extreme value 
statistics is indeed of the Gumbel type p5[| , while some corrections appear in the marginal case r] = 1. 
Interestingly, a related conjecture was proposed recently for the d = 2 problem with r] = 0, which 



corresponds to the localization of electrons in a random magnetic field |13|. In this two dimensional 
case, the choice rj = 2 — d = corresponds to a marginal logarithmic growth of the correlations. 

Returning to the one dimensional case, the situation changes drastically when r] < 1, which 
corresponds to a typically 'growing' profile E{x). The ratios c„ diverge with the system size for all 
temperatures, suggesting indeed a change of universality class. The only known possibility at present 
is then to describe the system within a 'continuous' rsb, which can be interpreted as a recursive tree- 
like construction of the low-lying energy state. In particular, the correlation of the low-lying states 
have a well known ultrametric structure. How well this ultrametric structure (known to be exact for 
the case where the dimension of x is infinite) reproduces the distribution and correlations of the low 
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lying states in finite dimension is an open problem Q We leave this problem for further studies [24|. 

In summary, we have argued for the one dimensional problem that the one step replica symmetry 
breaking scheme encodes exactly the results on the statistics of extremes for variables which are: 

i) not too correlated (i.e. when the above exponent rj is larger than 2 — d) 

ii) distributed asymptotically as generalized exponentials (i.e. as exp— 1£'|''). 

The case of long range correlations may correspond, in some limit, to the 'continuous' RSB scheme. 
However, the precise link between the two is not clear to us at present, and we think that it is an 
important path to explore further. It would also be interesting to think of the spin-glass problem from 
the point of view of the classification of very low energy states and excitations, which could perhaps 
provide a natural link between replicas and 'droplet-like' descriptions [26|. 



ACKNOWLEDGMENTS We thank A. Baldassarri, R. Cont and V. Dotsenko for many interesting 
discussions. MM thanks the SPhT in the CEA Saclay for its hospitality. 



Appendix 

In this appendix we explain briefly how the replica method and the direct probabilistic analysis 
lead to the same result for the two point correlation function in decaying Burgers turbulence at large 
times. We are interested in the case where E{x) has local correlations (see section 5.3), in which case 
the result of Kida reads: 

v{x)v{x + r) = vl,,H y^'j (82) 
where i{t) is the mean distance between two shocks, i^max = i{t)/t, and 

H[x) = —=—x / — (83) 

V27r dx Jo (p[x -I- y) -I- (p{x - y) 



and (j) is an error function: (j){x) = /q°° dzexp{—z + \J'k/2 xz). 

We shall sketch how these results can be obtained from the replica representation (|7^. The 
computations are lengthy and already contained in some previous papers. Here we just want to help 
the interested reader to find his way in the litter ature in order to obtain the result. One starts from 
the replicated partition function ([7^): 



1 " 

Z{xi)Z{x2)...Z{Xn) = YeiCl>[- ^ Rn{a),TT(h)XaXb] , (84) 



where the Rab matrix is a one step RSB matrix Q with elements Rab = f when a and b are in the same 
diagonal block of size m, Rab = when a and h are in different blocks, and Raa = (1 ~ m)r, enforcing 
the sum rule Ylb=i Rab = 0. The first step, derived in the appendix D of ||2^, deduces from ( |8^ ) the 
correlation between the powers n/2 of the partition function: 



Z(x,t)"/2Z(y,t)«/2 = „, I r^/.-"/^ / d.(e. + Me, 

B{—n/2, —n/2) Jo fi i 2tt J 



n/m 

■ (85) 



^Note that the average ground state energy predicted by the Gaussian variationnal replica theory does not lead back 
to the exact result |23] in the soluble random walk case rj — 0. 
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where we have defined 



-mr{z—y)'^/2 



Using the general link v = ~2i^ ^ , one gets 



v{x)v{y) = lim 



Z(x,t)"/22'(y,t)"/2 = {2vmrY{g 



11 + 912} 



n^o 'n? dxdy 

where we have defined (the notations are those of |20| | -Appendix B): 

/ dz {ex + /iCj/)™""^ (x - z)ex{y - z)ey 



911 = {l-m) dfi 



and: 



912 = m 



J dz {Cx + fJ-Cy)"^ 

(^J dz {ex + IJ-ey)""'^ {x - z)ex^ (^J dz (e^ + tJ.ey)"^~^ {y - z)e^ 



(86) 



(87) 



(88) 



(89) 



J dz {ex + f^ey)'^ 

This expression could also be derived directly without replicas from the infinite sum (|7^), with the 
identification m = fi, r = K/{Tfi) as can be seen from formulas (B10,B11) of (where the number 
l/(mr) was called 6). 

The whole problem is now to evaluate this expression in the limit of large Reynolds, which means 
small /X or low T. In this regime, using the fact that m scales linearly with T and r scales as l/T^, it 
has been shown in the appendix B of pO] that expression (p^) reduces to: 



v{x)v{y) = 




dh- 



,-hV2 ^^2 _ ^2/4] 



vr Jo 



dh- 



where 



Mo{x) 



e-M/2_y^o(/i_|) + e'^'^/27Wo( 



-z^/2 



2tt 



d 



\x - y\ 



(90) 



(91) 



So the natural length scale appearing in this solution is ^ = ^/T/Jjjlk). Using (73,39), one sees easily 
that it precisely scales at large times as the average distance between shocks of Kida's analysis. In 
terms of reduced lengths, one can check that the two distributions (pO|) and (p^) are actually identical. 
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